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\o ' 

I Aims. We present the results of multi epochs imaging observations of the companion to the planetary host Gliese 86. Associated to radial 
\^ velocity measurements, this study aimed at characterizing dynamically the orbital properties and the mass of this companion (here after Gliese 

86 B), but also at investigating the possible history of this particular system. 
\^ Methods. We used the adaptive optics instrument NACO at the ESO Very Large Telescope to obtain deep coronographic imaging in order to 

determine new photometric and astrometric measurements of Gliese 86 B. 

Results. Part of the orbit is resolved. The photometry of Gl 86 B indicates colors compatible with a ~ 70 Jupiter mass brown dwarf or a 
^ white dwarf. Both types of objects allow to fit the available, still limited astrometric data. Besides, if we attribute the long term radial velocity 
I residual drift observed for Gl 86 A to B, then the mass of the latter object is =^ 0.5 Mq. We analyse both astrometric and radial velocity data 
O to propose first orbital parameters for Gl 86 B. Assuming Gl 86 B is a ^ 0.5 Mq white dwarf, we explore the constraints induced by this 
hypothesis and refine the parameters of the system. 

• 

^ ' Key words. Stars individual: Gliese 86 - stars: low mass, brown dwarfs - planetary systems 

>y , 1 . Introduction mass) planet Gl 86 b, orbiting Gl 86 A at ~ O.UAU. This star 

is also surrounded by a more distant companion Gl 86 B, dis- 

One of the biggest challenges of today astronomy is to detect j <. on ait • u i j * j 

°° ° covered at ~ 20 AU using coronagraphy coupled to adaptive 

and characterize extra solar planetary systems, and to under- ... • /ci » i ^AA1^ t-u » j u » » t 

^ ■' ■' optics imaging (Els et al. 2001). The estimated photometry of 

stand the way(s) they form and evolve. Over the past decade, qj g . ^ compatible with that expected for a 40-70 Mj brown 

the technical improvements have allowed detections of more j ^ • o iv/r p m u - /onnc\ 

^ dwarf companion. Howewer, Mugrauer & Neuhauser (2005) 

than 150 extrasolar planets via radial velocity (hereafter RV) u j ..u- i . ui vu i u-. 

^ J ^ ' showed recently that this was also compatible with a cool white 

measurements down to 7.5 Earth Masses (Rivera et al. 2005) j f j 4.1. » i u »u • i-i 1 j 

dwarf, and that the latter hypothesis was more likely regard- 
around solar type stars, while direct imaging allows now the • 1^ u j . t • t-u u t 
■^^ ° ° ing the K band spectrum of the companion. The absence of 
detection of giant planets around young stars (Lagrange & m 1 1 j . • v n . u j 
° ^ ^ & \ b b near-IR molecular and atomic lines as well as the steep K-band 
Moutou 2004; Chauvin et al. 2004). From the theoretical point ^. ■ a a • . . -.u u . • ^ ac 

' ^ continuum are indeed consistent with what is expected for a 

of view the influence of the multiplicity or companionship with ,. , u - » -.u cc ^- ^ » i, - u »u 

^ -' ^ ^ high gravity object with an ettective temperature higher than 

outer bodies (e.g. brown dwarfs; hereafter BD) on the dynamics 4qqq 



and orbital stability of the inner planets has been highlighted. 
This has led to constant efforts trying to identify outer compan- 
ions to those stars hosting planets plus long term RV drifts. 

Gl 86 A is a KOV star with an estimated mass of 0.8 
(Siess et al. 1997; Baraffe et al. 1998) and is located at 10.9 pc 
from the Sun (Ferryman et al. 1997). Through RV measure- 
ments, Queloz et al. (2000) have detected a 4 Mj (minimum 



Apart for the RV wobble due to the hot Jupiter compan- 
ion, Gl 86 A also exhibits a long term RV drift measured with 
CoRAVEL and Coralie over 20 years. This drift indicates the 
possible presence of an additional more distant companion, 
with a substellar mass and a distance to star greater than ^ 20 
AUs. Els et al. (2001) claimed that Gl 86 B cannot account for 
this RV drift, due to its too low mass. They postulated instead 
Send offprint requests to: A.M. Lagrange that an additional companion, located in 2000 "behind" the star 

* Based on ESO observing programs 70.C-0543, 072.C-0624 and (i.e., under the coronographic mask), could be responsible for 
073.C-0468 at the VLT the observed drift. 
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Table 1. Observation log. NDsij^, is a CONICA neutral density filter with a transmission of 1.4%. S13 and S27 are two CONICA cameras 
corresponding respectively to a platescale of 13.25 and 27.01 mas. WFS corresponds to the wave front sensor of the adaptive optics system. 



UT Date 


Filter 


Camera 


Observation type 


Exp. Time (s) 


WFS 


Obs-Program 


Platescale calibrator 


12/11/2003 


Ks 


S27 


coronagraphy (0.7") 


100x0.6 


VIS 


072.C-0624 


01 


OriC 


12/11/2003 


2.17 + ND,,„rt 


S27 


direct 


15 X 4.0 


VIS 


072.C-0624 


01 


OriC 


22/09/2004 


H 


S13 


coronagraphy (0.7") 


48 X 1.0 


VIS 


073.C-0468 


01 


OriC 


22/09/2004 


H + NDsho,-, 


S13 


direct 


42 X 0.35 


VIS 


073.C-0468 


01 


OriC 


29/07/2005 


Ks 


S27 


coronagraphy (0.7") 


400 X 0.8 


VIS 


075.C-0813 


01 


OriC 


29/07/2005 




S27 


direct 


400 X 0.35 


VIS 


075.C-0813 


01 


OriC 


29/07/2005 


H 


S13 


coronagraphy (0.7") 


360 X 1. 


VIS 


075.C-0813 


01 


OriC 


29/07/2005 


H + ND,ho,-. 


S13 


direct 


400 X 0.35 


VIS 


075.C-0813 


01 


OriC 


29/07/2005 


J 


S13 


coronagraphy (0.7") 


165 X 2. 


VIS 


075.C-0813 


01 


OriC 


29/07/2005 


J + NDshort 


S13 


direct 


240 X 0.5 


VIS 


075.C-0813 


01 


OriC 




Fig. 1. VLT/NACO K^-band coronagraphic image of Gl 86 A and B, 
acquired on September 24, 2004, with an occulting mask of diameter 
0.7". 

In the course of a deep search for faint outer companions 
to stars hosting planets with NACO, we were able to make new 
images of Gl 86 A and B in the near IR. We present the obser- 
vational results in Sect. 2. In Sect. 3, we report new photometric 
result of Gl 86 B and we present an analysis of both astromet- 
ric and RV data, assuming that the RV drift is due to Gl 86 B. 
Finally, in Sect. 4 we discuss the nature of Gl 86 B, and we 
confirm that it is very probably a ~ 0.5 Mq white dwarf. We 
discuss the implications of this hypothesis. 

2. Observations 

2.1. NACO observing log 

Observations of Gl 86 were performed on November 12, 2003, 
September 22, 2004 and uly 29, 2005 with NACO at the VLT. 
NACO is equipped with an adaptive optics system (Rousset et 
al. 2002; Lagrange et al. 2002) that provides diffraction limited 
images in the near infrared (IR) and feeds the observing camera 
CONICA (Lenzen et al. 2002). Both coronagraphic and direct 
images were performed to image respectively Gl 86 B and A 



within the linearity domain of the detector. Note that between 
the two observing dates, the CONICA detector was changed 
and the latter detector was more efficient. 

The calibrations of platescale and detector orientation were 
done using the ©i Ori C astrometric field on November 12, 
2003, September 22, 2004 and July 29,2005. On November 12, 
2003 and July 29, 2005, the orientation of true north of the 
S27 camera was found respectively at -0.06° and -0.05° east 
of the vertical with an uncertainty of 0.20° and the platescale 
was 27.01 + O.lOmas. On September 22, 2004, the orientation 
of true north of the S13 camera was found 0.20° east of the 
vertical with an uncertainty of 0.20° and the platescale was 
13.25 + O.lOmas. Table 1 summarizes the new observations 
as well as archival ones and Fig. 1 shows a Kj (A^ = 2.2 yum, 

- 0.35 yum) image recorded in September 2004. 

2.2. Photometric measurements 

On July 29, 2005, the NACO measurements of G186 A and B 
were obtained in J,H and Ks filters, under photometric condi- 
tions. The JHKs contrasts were determined using the decon- 
volution algorithm of Veran & Rigaut (1998). Based on the 
JHK photometry of G186 A from the 2MASS All-Sky Catalog 
(Cutri et al. 2003), we then deduced the JHK photometry of 
G186B (see Table 1). The transformation between the K^ filter 
of NACO and the K filter used by CTIO-2MASS was found to 
be smaller than 0.03 magnitude. 



Table 2. Photometry of G186 A and B 



Component 


J 




H 




K 






(mag) 




(mag) 




(mag) 




G186A° 


4.79 ± 


0.03 


4.25 ± 


0.03 


4.13 ± 


0.03 


G186B* 


14.7 ± 


0.2 


14.4 ± 


0.2 


13.7 ± 


0.2 


G186 B' 


12.9 ± 


0.3 


13.1 ± 


0.2 


12.8 ± 


0.2 



" from the 2MASS All-Sky Catalog (Cutri et al. 2003). 
from Els et al. (2001) 

from " and NACO measurements presented in this work. 
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Fig. 2. Least square fits of the right ascension (left) and of the declination (right) of Gl 86 B relative to A. 



The reported photometry is significantly different from the 
ones given by Els et al. (2001) (see Table 1). This could be 
related to systematic photometric errors induced by an incor- 
rect subtraction of the Gl 86 A PSF wings within the coro- 
nagraphic image. This effect is generally more important for 
shorter wavelengths where AO corrections is poorer. It often 
leads to an underestimation of the companion flux, as this is 
the case in Table 1 when comparing ADONIS and NACO data. 
Thanks to the highest angular resolution and the enhanced de- 
tection capabilities provided by NACO at VLT, we can reason- 
ably expect to be less sensitive to this PSF subtraction effect to 
derive the Gl 86 B photometry. 

The new NACO photometry found is still compatible 
with the conclusions of Els et al. (2001) that G186B has a 
photometry similar to that expected for a substellar companion 
with a mass of 40-70 Mj (spectral type L7-T5). However, 
this photometry can also correspond to the one expected for 
a cool white dwarf and, as recently claimed by Mugrauer & 
Neuhauser (2005), this is more likely the case as the spectrum 
of Gl 86 B does not exhibit the molecular absorption features 
in K band that are characteristic for L or T dwarfs. 

In the following, we reinvestigate this issue (brown or white 
dwarf) from a dynamical point of view. 



2.3. Astrometric measurements 

The offset positions of Gl 86 B to A, recorded with NACO 
on 12 November 2003, 22 September 2004 and 29 July 2005, 
were translated into physical values using the corresponding 
astrometric calibration data. The shifts induced by the use of 
different filters between coronographic and direct images were 
taken into account. Table 3 suimnarizes the measured values 
and Figure 2 shows the various data points in a (Ac, A6) dia- 
gram, as well as the offset positions of Gl 86 B to A measured 
by (Els et al. 2001) with ADONIS/SHARPH on 8 September 
2000. The orbital motion of Gl 86 B is clearly identified. This 
confirms the independent detection of Mugrauer & Neuhauser 
(2005). 



Table 3. Offset positions of the Gl 86 B relative to A 



UT Date Juhan Aor 


A<5 


Separation Position 


Date (mas) 


(mas) 


(mas) Angle (°) 



08/09/2000 2451796 1510 + 25 -853 + 6 1734 + 22 119.5 + 0.8 

10/11/2000 2451859 1522 + 3 -789 + 20 1714 + 10 117.4 + 0.4 

12/12/2000 2551891 1508 + 10 -851 + 13 1732+ 11 1194 + 0.4 

12/11/2003 2452986 1818 + 12 -573 + 4 1906 + 11 107.5 + 0.4 

22/09/2004 2453271 1872 + 16 -513 + 6 1941 + 14 105.3 + 0.5 

29/07/2005 2453581 1920 + 13 -435+9 1969 + 11 102.7 + 0.4 



2.4. Radial velocity data 

Radial velocity measurements of Gl 86 A have been gathered 
for more than 20 years now. The whole data set reveals in ad- 
dition to a short period modulation of ~ Ikms"^ ampUtude 
that has been attributed to a hot Jupiter companion (Queloz 
et al. 2000), the presence of a regular continuous decrease of 
~ 2kms~^ in 25 years (Fig. 7). 

It is tempting to try to attribute this regular decrease to 
Gl 86 B. The temporal derivative of the radial velocity of the 
primary in a binary system is easy to derive. One gets 



dVr 



Gm 



—— = T- sin i sin(w + v) 

at 



(1) 



where G is the gravitational constant, m is the mass of the com- 
panion r is the distance between the two bodies, i is the incli- 
nation of the orbit with respect to the plane of the sky, a> is 
the argument of periastron, and v is the true anomaly, i.e. the 
current polar position along the orbit with respect to the perias- 
tron. Of course most of these quantities are unknown, but a sim- 
ple application assuming sin / sin(w + v) ^ 0.5 and r ^ 20 AU 
shows that dvr/df -~ 2kms V25 yr is hardly compatible with 
m = 70 Mj, but rather with m ranging between 0.2 and 1 Mq. 

This result led Els et al. (2001) to conclude fliat the RV 
residuals are not due to Gl 86 B, but rather to an unseen, addi- 
tional body. Conversely if we keep attributing the RV decrease 
to Gl 86 B, this raises the question of the mass of Gl 86 B. The 
available photometry is compatible with a 70 Mj object (Els et 
al. 2001). But it can also be compatible with a ~ 0.5 M© object 
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if this object is a white dwarf. Obviously, more data, in particu- 
lar spectroscopic data are needed to discriminate between these 
two possibiUties. 



3. Data Analysis 

3. 1. General analysis of astrometric data 

From Fig. 2, one can see that on the plane of the sky, the four 
points (see plots below) are roughly aligned, so that the only 
relevant information we can derive from these data is a middle 
astrometric position (at t = 2003.0) and temporal derivatives 
of the right ascension a and of the declination 6. We thus per- 
form a least- square fit of the available data to derive them. The 
result is shown on Fig. 2. We see that a and 6 actually vary 
roughly linearly with time. The linear fit is therefore relevant. 
The corresponding temporal derivatives are 



10 



— = -89.5 ±8.7 mas yr"' 
df 

m 

dt 



85.6 ±7.18 mas yr"^ 



(2) 



These derivative values, together with the mean present values 
of a and 6, provide four contraints on the orbit of the compan- 
ion with respect to the primary. In principle, this orbit is fully 
characterized by 6 orbital elements, plus the unknown mass m 
of the companion. The constraints allow us to fix 4 of them. We 
chose to let the mass m of the companion, the inclination / with 
respect to the plane of the sky, and the longitude of the ascend- 
ing node Q as free parameters. For any given set of parameters 
(m, i, O) we are able to derive the remaining ones, i.e. the semi- 
major axis a, the eccentricity e, the argument of periastron w 
and the mean anomaly M. We recall that M is a quantity that 
characterizes the present position of the companion on its orbit. 
M is proportional to the time, M = at periastron and M = 2n 
one orbital period later. 

3.2. Analysis assuming that Gl 86 B is a 70 Mj object 

Depending on the free parameter set we choose, there is not 
necessarily an orbital solution compatible with the constraints. 
In particular, it turns out that there is no solution for i < 120°. 
This means that we are viewing the orbit nearly from its south 
pole. The result of the parameter space exploration is shown in 
Figs. 3-5. The semi-major axis, the eccentricity, and the mean 
anomaly are plotted as a function of O, for different values of 
the inclination /, and for a fixed companion mass m = 70 Mj. 
We note that in some cases (i = 120°) there is not a solution 
for every Q value. We note also that the orbit is necessarily 
eccentric (e > 0.35 in any case), and that in all cases, the com- 
panion is presently short after periastron (0 < m < 60°). Of 
course we explored other companion masses in the compatible 
range (60 Mj ^ m ^ 90 Mj). The result is not shown here but it 
is nearly equivalent to that for m - 70 Mj. Actually Figs. 3-5 
represent the standard solution. 
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Fig. 3. The semi-major axis a of the orbital solution for the Gl 86 
companion, as a function of the longitude of the ascending node CI, 
for various values of the inclination i between 120° and 180°, for a 
fixed companion mass m = 70 Mj. The bullet represents the solution 
plotted in Fig. 6 (upper plot) and detailed in Eq. (3). 
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Fig. 4. Same as Fig. 3, but for the orbital eccentricty of the solution 



In order to better show the shape of the orbital solution, we 
display one particular, typical solution, marked as a bullet in 
Figs. 3-5, and charaterized by ; = 150°, and 



a = 47.58 AU , e = 0.6185 , Q. = 300° , 
w= 19.71°, M= 18.58° 



(3) 



The projection of this solution onto the plane of the sky is 
shown in Fig. 6. We clearly see that the orbit is eccentric and 
that the present day position of the companion is short after pe- 
riastron. The associated orbital period is 353 yr, and the last pe- 
riastron passage turns out to have occurred in 1984. Of course 
the latter quantities are subject to some variations if we con- 
sider another solution. 

In Fig. 7, we show the Gl 86 radial velocity data set, su- 
perimposed to the theoretical curve that would be expected for 
the solutions we display in Fig. 6. Note that in those curves, we 
dot not add the short period modulation due to the hot Jupiter 
companion, as this object produces a much smaller amplitude. 
We also add to the theoretical radial velocity curve an empirical 
offset, intended to correspond to the mean heUocentric velocity 
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Fig. 6. A representation of the orbital solutions described by Eq. (3) (upper plot, brown dwarf case) and Eq. (4) (lower plot, white dwarf case), 
as projected onto the plane of the sky. In each case. The right plot represents a view of the full orbit and the left plot is an enlargement of the 
present day motion. The dotted line is the projection of the line of apsides of the orbit. 



of the Gl 86 system, fixed in such a way that the radial velocity 
matches the mean observed value in 2003.0 . Actually the only 
relevant parameter we need to compare between the data and 
the theory is the mean temporal derivative of the radial velocity 
at in 2003.0, and also the general trend over 25 years. 



In Fig. 7, the theoretical radial velocity curve correspond- 
ing to Eq. (3) is represented as a dashed line. We see that it does 
not match the data. In fact the decrease in 2003.0 is only 10% of 
the observed values (0.1 kms"V25 years). As explained above, 
this was expected from our order of magnitude estimate of the 
mass needed to account for the observe decrease rate. 



3.3. Analysis assuming Gl 86 B isa~ 0.5 Mq object 



If we now assume that the residuals of the radial velocity data 
are due to Gl 86 B, we get additional constraints to the orbital 
parameters. In particular, we can force the temporal derivative 
of the radial velocity in 2003.0 to match the observed one. This 
in turn enables to fix the mass m of the companion instead of 
giving it as an input parameter. However, this single criterion 
turned out not to be sufficient. We may derive solutions that 
fit the radial velocity derivative in 2003.0 but that do not fit 
the radial velocity data over the whole observation period, es- 
pecially the older data. Hence we retain in the fitted solutions 
only those which fit a convenient least square criterion with the 
whole radial velocity data sample. 
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Fig. 7. Radial velocity data of Gl 86 A as monitored over 25 years, 
superimposed to the theoretical curves (in grey) corresponding to 
the orbital solutions displayed in Fig. 6. The low accuracy data up 
to 1998 are the Coravel data (typical error ±0.27 km s"'), while the 
subsequent high accuracy data are the Coralie ones (typical error 
+0.005 kms"'). The dashed curve corresponds to the orbit detailed 
in Eq. (3) where Gl 86 B is taken as a brown dwarf. The solid curve 
corresponds to the orbit described in Eq. (4), where Gl 86 B is fitted 
as a white dwarf. The fit of the radial velocity residuals is much better. 



The result of the exploration of the parameter space is 
shown in Figs. 8-11. Note that contrary to Figs. 3-5, solutions 
are plotted only for -83° < O < 10°; actually there is no con- 
venient solution out of this range of Q. We see also that there 
are solutions now for 1 10° < i < 150°. The orbit is still viewed 
from south but it does not exactly lie in the plane of the sky 
{i = 180°). Actually with exactly i = 180°, there would be 
no radial velocity signature. The significant decrease of the ra- 
dial velocity as observed over 25 years forces the inclination / 
not to be too close to ; = 180°. The solutions are still eccen- 
tric, and the present location of Gl 86 B is still more or less 
short after periastron. The most interesting outcome concerns 
the now fitted mass of the companion (Fig. 11). No solution 
with m < 0.4 Mq is found, and the more likely solutions corre- 
spond to 0.4 Mq < m < 0.6 Mq. This is of course very different 
from typical brown dwarf values, but falls well in the range of 
typical white dwarf masses. 
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Fig. 8. The semi-major axis a of the orbital solution for Gl 86 B that 

also fits the radial velocity data residuals, as a function of the longi- 
tude of the ascending node CI, for various values of the inclination / 
between 110° and 150°. The bullet represents the solution plotted in 
Fig. 6 (lower plot), described in Eq. (4). 
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Fig. 9. Same as Fig. 8, but for the orbital eccentricty of the solution 



As in the previous section, we display one peculiar solution 
assumed to represent a standard solution, charaterized by i = 
150° and 



a=18.42AU, e = 0.3974 , Q = -35° , 
0) = 18.05° , M = 100.5° , m = 0.4849 Mq 



(4) 



This solution is marked as a bullet in Figs. 8-11. The orbital pe- 
riod now only 69.7 yr, and the last periastron passage occurred 
in 1983. 

In Fig. 6, we show the projection of this solution onto the 
plane of the sky like we did it for the orbit corresponding to 
Eq. (3), and in Fig. 7 we show the correspondant radial velocity 
curve as a solid grey curve. The agreement with both the radial 
velocity and the astrometric data is very good. Apart from small 
changes in the orbital elements, the main difference with the 
orbit described in Eq. (3) is the mass of the companion. With 
m = 0.5 Mq, it is obviously not a brown dwarf. 
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Fig. 11. Same as Fig. 8, but for tlie fitted mass of the Gl 86 B compan- 
ion 



4. Discussion 

4.1. The nature of Gl 86 B 

From the above analysis, either Gl 86 B is a brown dwarf, and 
then it is unable to explain the RV residuals, either it is a ~ 
0.5 Mq white dwarf. In the former case, another, massive object 
is required to explain the RV residuals. In that case, one should 
wonder why this object has not been detected yet, unless it is 
angularly close to the primary, so that it should disappear under 
the coronographic mask used in the images, as suggested by Els 
et al. (2001). Given the inclination we derive for Gl 86 B, the 
whole system is thus far from being planar. Independently from 
the low probability that such an additional massive component 
would be located today in such a position that it could not be 
detected, the dynamical stability of the whole system should 
be questioned. It is well known (Beust et al. 1997; Beust 2003; 
Krymolowski & Mazeh 1999) that multiple systems with high 
mutual inclinations are often subject to the Kozai resonance, 
and that this can lead to instability. 

It seems thus more natural to try to attribute the RV resid- 
uals to the sole Gl 86 B companion. In that case, it must be a 
0.4-0.6 Mq object. As from its photometry it cannot be a main 



sequence star of that mass, Gl 86 B turns out to be necessar- 
ily a white dwarf. Our dynamical analysis finally leads to the 
same conclusion Mugrauer & Neuhauser (2005) derived from 
independent spectrophotometric arguments. 

Based on the present constraint put on the mass of Gl 86 B 
and on the new NACO JHKs photometry, presented in Sect. 
2.2, we can now re-investigate the physical properties of this 
white dwarf companion, using predictions of the evolution- 
ary cooling sequences models of Bergeron et al. (2001) for 
hydrogen- and heUum-rich white dwarfs. 

The model predictions are reported in a color-magnitude 
diagram (J - K vs M^r) for both cases: hydrogen-rich (Fig. 12, 
left) and helium-rich (Fig. 12, right) white dwarfs. We can first 
notice the discrepancy between the model predictions and the 
previous photometric data of Els et al. (2001) that Mugrauer 
& Neuhauser (2005) used to derive an effective temperature of 
5000 ±500 K for Gl 86 B. Our new NACO photometric data are 
in very good agreement with the model and with the dynamical 
constraints. Then, if we add the fact that the mass of Gl 86 B is 
dyanamically constrained between (0.4-0.6 Mq), we can derive 
the effective temperature, the gravity as well as the cooling age 
of the Gl 86 B companion based on models predictions. The de- 
rived physical parameters for hydrogen- and helium-rich white 
dwarf model predictions are reported in Table 4. 

4.2. The Initial-Final Mass Relationship 

Both dynamical and spectrophotometric studies come come to 
the conclusion that Gl 86 B is actually a white dwarf. Let us 
now investigate the dynamical implications of this hypothe- 
sis. The main uncertainty concern the initial main-sequence 
mass of Gl 86 B before its evolution to the white dwarf 
state. This general problem is known as the Initial-Final Mass 
Relationship (IFMR) for white dwarfs (Jeffries 1997). This 
problem, together with the upper mass limit for white dwarf 
progenitors, has been the subject of intense investigations in 
the past (Weidemann 1977, 1987, 1990). 

For what concerns Gl 86, a first constraint is that in any 
case, Gl 86 B must have been more massive than Gl 86 A in the 
past (i.e., 0.8 Mq), in order to have more quickly evolved to the 
post main-sequence state. 



Table 4. Physical parameters of Gl 86 B based on predictions of the 
evolutionary cooling sequences models of Bergeron et al. (2001) for 
hydrogen- and helium-rich white dwarfs. 



Model 


Mass 
(Mo) 


Teff 

(K) 


log(g) 


Cooling Age 
(Gyr) 


H-rich 


0.4 
0.5 


5500 ± 1000 
6000 ± 1000 


7.66 + 0.02 
7.86 + 0.01 


1 4+1'* 

^- -0.42 

^Kl 




0.6 


7000 ± 1000 


8.01 ±0.01 




He-rich 


0.4 


6000 ± 1000 


7.70 ± 0.01 


1 6*^-^ 

■^•"-0.6 




0.5 
0.6 


7000 ± 1000 
8000 ± 1000 


7.88 ± 0.01 
8.03 ± 0.01 


■■'-0.34 
1 9+0.6 
• -0.29 
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-10 12 -10 12 

J-K J-K 



Fig. 12. Color-magnitude diagram (J-K vs M^:) with model predictions for different masses in two different cases: white dwarfs with hydrogen- 
rich (Fig. 12, left) and helium-rich (Fig. 12, right) atmospheres. The predictions for a 0.4, 0.5 and 0.6 M© white dwarf, which is likely the case 
for Gl 86 B based on our dynamical analysis, are given in solid lines, the others in dotted lines. The iso-re(f lines {red) have been also reported 
in the zoom-in images. The photometric data from Els et al. (2001) and this work have been reported in both figures for direct comparison with 
models predictions. 



The IFMR is an increasing function of the initial mass. It 
is usually measured using white dwarfs that are members of 
open clusters of known ages. Weidemann (1987) gives a semi- 
empirical IFMR, but further measurements of white dwarfs 
in NGC 2516 (Jeffries 1997) have shown it was inaccurate. 
More relevant relations for various metallicities (Z) are given 
by Hurley et al. (2000). In the following, we will assume the 
IFMR given by Hurley et al. (2000) (Fig. 18) for Z = 0.02. 

Note that this IFMR is different from another one that is 
sometimes shown (Iben 1991; Bressan et al. 1993; Fagotto et 
al. 1994), which shows the mass of the white dwarf remnant as 
a function of that of the core at the beginning of the TP-AGB 
phase. We are interested into the full initial mass of Gl 86 B 
at Zero Age Main Sequence (ZAMS), so that the first IFMR is 
relevant here. 

4.3. Mass loss in a binary system 

Additional constraints can be derived if we consider the past 
evolution of the mutual orbit of Gl 86 A and B. The impor- 
tant post main-sequence mass loss of Gl 86 B that led it to 
its white dwarf state induced an evolution of the orbit that can 
be estimated. The general problem of orbital evolution due to 
mass loss in a binary system has been theoretically investigated 
by many authors. Basically, one must distinguish between slow 
and rapid mass loss. In the former case, the semi-major axis ap- 
pears to grow during the mass loss process, while the eccentric- 
ity remains secularly unchanged (Jeans 1928; Hadjidemetriou 
1963; Verhulst 1972); in the latter case (rapid mass loss) both 
the semi-major axis and the eccentricity grow (Blaauw 1961; 
Hut & Verhulst 1981). A major difference is that in the case of 
slow mass loss, the orbit always remains bound (il just widens). 



while in the latter case in can be disrupted. This actually occurs 
if the mass loss overcomes half of the mass of the whole system 
(Blaauw 1961). This case corresponds typically to supernovae. 

In the case of Gl 86, we shall be concerned by the slow 
mass loss case. The equations defining the variation of the 
semi-major axis a and of the eccentricity e are given by 
Hadjidemetriou (1963) : 

de M 

— - -{e + cos f) — ; (5) 
df ^' M ^ ' 

aGM{\ - e^) - constant , (6) 

where M is the total mass of the system, M the mass loss rate 
(due here to Gl 86 B only) and / is the true anomaly along the 
orbit. The second equation arises from the fact that the specific 
angular momentum C = r A v is unchanged. The first one is 
derived assuming that the change of the specific orbital energy 
U is only due to the the mass loss (df//df - -GM/r) where r 
is the radius vector (Verhulst 1974). 

If the mass loss is a slow process, Eq. (5) can be averaged 
over one orbital period. This turns out to give de/dt - 0, which 
means that the eccentricity is secularly constant (Jeans 1928; 
Hadjidemetriou 1963). Subsequently, the evolution of the semi- 
major axis just obeys the simple rule aM - constant. As M 
decreases, it is obvious that the orbit gets wider If the total 
change of M (only due to Gl 86 B) is known from the IFMR, it 
is then possible to derive the initial semi-major axis. 

4.4. Application to Gl 86 A and B 

If we apply this theory to the case of Gl 86 B, we are able to 
derive the former characteristics of the Gl 86 system. The fit of 
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Fig. 13. The initial ZAMS semi-major axis flinit of the Gl 86 B as com- 
puted for all solutions displayed in Figs. 8-11, using the IFMR from 
Hurley et al. (2000) and Ma = constant . Only the solutions that fit all 
the constraints have been retained (see text). 
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Fig. 14. Same as Fig. 13, but for the initial main-sequence mass of the 
Gl 86 B progenitor, assuming Gl 86 B is presently a white dwarf 
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Fig. 15. Same plot as Fig. 11, but all solutions laeding to unphysical or 
unacceptable values for a^^n e^^i, or nii^n have been removed. The grey 
shaded area corresponds to all possible values if we let the inclination 
i vary. 



Sect. 3.3 allows to derive the present day orbital and mass char- 
acteristics of Gl 86 B (a, e and m). For each solution, using the 
IFMR of Hurley et al. (2000), we are able to derive the initial 
mass mimt, and subsequently the initial initial semi-major axis 
flinit of the orbit, using aM = constant . All solutions that lead 
to unreasUstic (negative) values for Umit are then eUminated; 
we also eliminate all solutions for which mink < 0.8 Mq, as 
Gl 86 B must have been initially more massive than Gl 86 A. 
This can be done for every solution that fits the radial velocity 
and the astrometric data. This constraint turns out to be by far 
the strongest one. 

The result is shown in Figs. 13-14. In these figures, we plot 
the resulting values of aian, and numt for all the solutions dis- 
played in Figs. 8-11. However, we only retain those solutions 
which lead to compatible values for ainit, and to > 0.8 Mg. 
This is the reason why the curves are often interrupted. In par- 
ticular, all solutions with / = 110° have been eliminated. 

In all cases we have ainit < a (typically ai^t - 0.5a), show- 
ing that the orbit is more detached presently than it was in the 
past. This is for instance the case for the solution described in 
Eq. (4), for which we have 



= 12.97 AU 



= 0.865 Mp 



(7) 



This solution is marked as bullets in Figs. 13-14. We see 
that ainit is not very strongly constrained. The original mass 
of Gl 86 B is better constrained. On Fig. 14, we see that it 
may range between 0.8 and 2 M©, but more probably it was 
< 1.5 Mq. In fact, the solution giving minit - 2Mo are those 
which correspond to the smallest values for Ainit (see Figs. 13- 
14). If ainit was too small, the past orbital stability of the ex- 
oplanet companion of Gl 86 A may be questioned. Obviously 
this dynamical issue needs to be investigated into further de- 
tails. But as a first attempt, let us consider a possible origi- 
nal configuration of Gl 86 with a 0.8 Mq Gl 86 A and a 2 Mq 
Gl 86 B progenitor. The Hill radius around Gl 86 A can thus 
be estimated to ~ 0.45 d, if d is the separation between the two 
stars. If we take for d the periastron of the orbit, with e ^ 0.3 
(this is the value derived for such solutions; see Fig. 9), and if 
we assume that the Hill radius must be at least ~ 2 times larger 
than the 0.11 AU semi-major axis of the planet to ensure sta- 
bility, we derive ainit ~ 0.7 AU; actually for all solutions with 
ainit < 1 AU, the orbital stability of the exoplanet is subject to 
caution. 

Another puzzling issue is the way the exoplanet formed. To 
what extent was the initial circumstellar disk of Gl 86 A that 
gave birth to its companion truncated by tidal interaction with 
Gl 86 B ? According to Eggenberger et al. (2004), the mini- 
mum separation in a binary that allows a large enough circum- 
stellar disk for planet formation to survive ranges between 10 
and 50 AU. This could mean that we should remove all solu- 
tions with ainit < 10 AU, which would result in minit < 1-3 Mq. 

The constraints on ainit and minit help to eliminate some of 
the fitted solutions in Figs. 8-11. This does not change the basic 
constraints on a, e and M, but refines that on the present mass 
m of Gl 86 B. In Fig. 15, we show the same plot as in Fig. 11, 
but all solutions that do not fulfill the constraints on ainit, ^init 
and minit have been removed. In order to explore all possilities, 
we performed the same calculation for many inclination values 



10 



A.-M. Lagrange et al.: New constrains on GUese 86 B 



(not only for i = 120°, i = 130°. . . ). The resulting possibilities 

are summarized as grey areas in Fig. 15. We see that m is fairly 
well constrained. It is thus possible to definitively stress that 

0.48 Mo < m < 0.62 Mo (8) 

and even probably we could say that m < 0.55 Mg. The sharp 
lower limit at m = 0.48 M© is due to the lower limit of 0.8 Mo 
for OTinit; the upper Umit at m — 0.61 Mq corresponds to Uiiat — 
0. 

5. Conclusion 

The identification of the orbital motion of Gl 86 B around 
Gl 86 A, combined to the measured residuals of the radial ve- 
locity data, allow to severely constrain the whole Gl 86 sys- 
tem and its past evolution. Our dynamical study shows that 
Gl 86 B is very probably a white dwarf, in agreement with the 
conclusions of totally independent spectrophotometric study by 
Mugrauer & Neuhauser (2005). The brown dwarf hypothesis of 
Els et al. (2001) can therefore be definitively ruled out. 

The mass of Gl 86 B is severely constrained by the dynam- 
ics. We derive 0.48 Mq < m < 0.62 Mq. The orbit is eccentric 
(e > 0.4) with a semi-major axis of a few tens of AU. The 
associated orbital period is several hundreds of years at least, 
and the stars have recently (5-20 years ago) passed at perias- 
tron. The orbit is retrograde with respect to the plane of the sky, 
but does not exactly lie in that plane. Actually we can say that 
120° <, i ^ 150°. 

Based on new photometric results on Gl 86 B and the dy- 
namical mass constrains, we also re-investigated the physical 
properties of this white dwarf companion. Using model predic- 
tions of Bergeron et al. (2001), we derived the effective tem- 
perature, the gravity and the cooling age of Gl 86 B for both 
hydrogen-rich and helium-rich atmospheres models of white 
dwarfs. 

When Gl 86 B was a main sequence star, its mass probably 
ranged between 0.8 M© and 1.5 M©, which impUes a spectral 
type between K2V and F7V. Its orbit was closer. The strong 
post-main sequence mass loss caused the orbit to widen. If it 
had been a more massive star, the initial semi-major axis would 
have been too small to allow orbital stability for the exoplanet 
orbiting Gl 86 A. 

However SafFe et al. (2005) recently used chromospheric 
index and metallicity measurements to estimate the age of all 
known stars harbouring exoplanets. For Gl 86 A, they derived 
an age ranging between 2 Gyr and 3 Gyr. Given the main se- 
quence lifetimes and the white dwarf cooling times (Table 4), 
assuming this age for Gl 86 B would imply that its progeni- 
tor had minit ^ 2 Mq. This seems to be incompatible with our 
dynamical constraints. Obviously, in order to solve this dis- 
crepancy, the dynamical evolution of the whole system, includ- 
ing the exoplanet needs to be investigated into further details. 
There are many open questions associated with this issue: the 
exoplanet must have survived all the late evolution stages of 
Gl 86 B. If the system is not coplanar, the exoplanet could have 
been subject to the Kozai resonance in the past. Moreover, the 
planet must have formed in a large enough circumstellar disk, 
which impUes a minimum initial separation of ~ 10 AU. All 



these issues need to be addressed, and this wiU be the purpose 
of forthcoming work. 
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